Multiplicity fluctuations in high energy hadronic and nuclear collisions 
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The showers of cosmic rays entering the Earth's atmosphere are main sources of information on cosmic rays and 
are also believed to provide information on elementary interactions at energies not accessible to accelerators. In 
this context we would like first to remind the role of inelasticity K and elementary cross section a and then argue 
that similar in importance are fluctuations of different observables. The later will be illustrated by multiplicity 
fluctuations in hadronic and nuclear collisions. 



1. INTRODUCTION 

The EAS (Extensive Air Showers) initiated 
by cosmic rays (CR) entering the Earth's atmo- 
sphere are our main tool to investigate the CR as 
such and they are believed to remain our only 
possible window to look at the elementary in- 
teractions at energies not accessible to acceler- 
ators 1 . Agreeing fully with the first part of 
this statement we would like to remind here that 
EAS are themselves macroscopic stochastic pro- 
cesses of variable length and this fact sometimes 
leads to results, which are unexpected from the 
elementary particle interactions point of view, the 
apparent intermittency seen there and connected 
solely with stochasticity of EAS only |2j being the 
best example. Other example worth to mention 
at this point is the apparent self-organized char- 
acter of EAS depending in a visible way on their 
length and clearly seen in data 0. Both phe- 
nomena do not depend on details of elementary 
interactions used in describing EAS. Actually, as 
it was widely appreciated some time ago |3J and 
is being rediscovered also at present out of 
numerous parameters of models of elementary in- 
teractions used in MC codes describing produc- 
tion and development of EAS the most important 
are: inelasticity K and elementary cross section 
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a (actually, the hadron - air nucleus cross sec- 
tion). Out of these two, the inelasticity K, i.e., 
the fraction of energy used for production of sec- 
ondaries and therefore not available for the sub- 
sequent interactions in developing EAS, is nowa- 
days not really a parameter as it used to be before 
(cf. 0]) but rather a number which must be re- 
produced by combination of all parameters used 
in the present models [S] 5 . 

2. FLUCTUATIONS 

However, there are some data (cf. [7]) which 
seem to demand more information. In particu- 
lar, in [H] we have shown that the phenomenon of 
the so called long- flying- component can be most 
naturally explained by allowing for the fluctua- 
tions in the cross section (which at that time 
were also widely discussed in the literature in 
other context). This conjecture was later quan- 
tify by using the so called nonextensive statisti- 
cal approach (essentially corresponding to using 
Tsallis entropy characterized by parameter q such 
that for q — > 1 it becomes the usual Boltzmann- 
Gibbs entropy) 9 . In such approach parameter q 
measures, in a sense, the amount of fluctuations. 

5 Actually K appears usually together with a, therefore, 
as we have shown in one needs at least two indepen- 
dent different experiments (measuring quantities in which 
K enters in different way) to be able to estimate K(s) 
from experiment in a model independent way. For other 
possibilities see |S| and discussion below. 
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For the purpose of this presentation it is enough 
to say that whenever in the exponential formula, 
exp(— X/A), parameter 1/A fluctuates according 
to gamma distribution then one should use in- 
stead expression 6 



exp (-X/A ) = [1 - (1 - q)X/A} 1 ^ 



(1) 



where 1/A = (1/A) and q = ((1/A) 2 ^ / (1/A) 2 

(averages are over the above gamma distribu- 
tion). We first applied this approach to the above 
mentioned data on long flying component an d 
found that fluctuations they correspond to are 
given by q = 1.3. 
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Figure 1. Results on the energy dependence of 
inelasticity K from pp and pp data (see |13| for 
details). 
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Figure 2. Results on the energy dependence of 
fluctuations as given by q — q(s) (points) com- 
pared with the best fit to energy dependence of 
parameter 1/fc of Negative Binomial multiplicity 
distribution (line, see text and ^Sj for details). 



6 See 1101 for details. Generalization of this approach to 
other type of fluctuations is discussed in 1111 . 



In Figs. 1-3 we provide results on inelastic- 
ity K and fluctuations seen in the rapidity dis- 
tributions, which originate from the multiplicity 
fluctuations |13I14| . They were obtained using 
nonextensive version of the information theory 
approach to multiparticle production processes 
based on the Tsallis g-entropy. Given the mean 
multiplicity (n) and mean transverse momentum 
(pt) one gets the most probably and least bi- 
ased rapidity distribution, being of the exp ? (...) 
type, i.e., depending on the nonextensivity pa- 
rameter q, which further depends only on the 
amount of energy available for production of sec- 
ondaries, W — K ■ y/s, i.e., on the inelasticity K. 
In this way inelasticity K and nonextensivity pa- 
rameter q specifying entropy used as measure of 
information are the only parameters to be fitted. 
Fig. 1 shows the energy dependence of inelastic- 
ity, K = K(s), whereas Fig. 2 the energy depen- 
dence of fluctuations as given by q = q(s). In the 
case when one knows rapidity distributions also 
for a given multiplicity (or, at least, for some mul- 
tiplicity bins) this method allows to deduce also 
inelasticity distribution, as can be seen in Fig. 
3 where it was done for ^/s = 200 GeV (black 
symbols, K = 0.52) and 900 GeV (open symbols, 
K = 0.38). They can be fitted either by gaussians 
(full lines) or, better, by lorentzian (dotted lines) 
shapes. 
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Figure 3. Examples of inelasticity distributions 
x{K) obtained using the available UA5 data on 
rapidity distributions in multiplicity bins (see ^21 
for details). 
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In what follows we shall, however, concentrate 
on fluctuations described by the parameter q 
shown in Fig. 2. We argue that these fluctua- 
tions are due to fluctuations of mean multiplic- 
ity used as our input when deducing inelasticity 
K. Namely, when there are only statistical fluc- 
tuations in the hadronizing system one expects 
Poissonian form of the corresponding multiplic- 
ity distributions. The existence of intrinsic (dy- 
namical) fluctuations means then that one allows 
mean multiplicity n to fluctuate. In the case when 
these fluctuations are given by gamma distribu- 
tion with normalized variance D(n) then, as a 
result, one obtains the Negative Binomial mul- 
tiplicity distribution, which depends on two pa- 
rameters: the mean multiplicity (n) and the pa- 
rameter k (k > 1) affecting its width, 



That is because (see [T5]1: 

P{n) = I dn : — — = 

Jo nl r W 

T(k + n) 7 fc 

r(i + n)r(fc) ( 7 + i) fe +" 1 ' 

where 7 = -Mr. According to our philosophy it 
is therefore natural to describe these fluctuations 
by the nonextensivity parameter q assuming that 
D(h) = q — 1, i.e., that q = 1 + 1/fc. However, in 
the nonextensive approach of this type q is lim- 
ited by q < 3/2 [T3|. This condition would then 
impose constraint on the amount of fluctuation 
as it corresponds to k < 2, the saturation limit, 
which from the naive extrapolation of tendency 
presented in Fig. 2 would appear at energies ~ 33 
TeV (or E LAB = 0.5 • I0 18 eV) range. Notice that 
this energy range of the ultra high energetic cos- 
mic rays in which effects connected with the GZK 
cut-off starts to be important . It means there- 
fore that the knowledge of fluctuations will most 
probably turn out to be as important as that of 
inelasticity and cross section. 

We would like to add to this one more observa- 
tion. In Fig. 4 we have plotted experimental data 
show on mean charged multiplicity (n c h) at differ- 
ent energies versus corresponding total inelastic 



cross section a (Fig. 4a) and the estimated vol- 
ume of the interaction region V (taken as being 
proportional to <r 3 / 2 ). The linearity of {n c h) and 
a is remarkable. This suggest that fluctuations of 
the multiplicity distributions should also be con- 
nected with fluctuations of the cross section but 
at the moment there are no investigations in this 
direction 7 . 
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Figure 4. (a) Experimentally observed correla- 
tion between (n c h) and cr. (b) The same but with 
a replaced by the volume of interaction estimated 
as being cx cr 3 / 2 . 

3. SUMMARY 

To summarize: phenomenon of EAS is macro- 
scopic and so complicated that it depends only 
on a small number of characteristics of elemen- 
tary interactions. It seems that they are: inelas- 
ticity K, cross section a and its fluctuations. We 
have argued that although we do not measure di- 
rectly the later we measure fluctuations of mul- 

7 Although constraints on the production cross section im- 
posed by the observed behaviour of multiplicity distribu- 
tions have been investigated already long time ago I17| . 
this subject is not pursued at moment. In models so far 
there is problem with (n c h) at energies of interest to CR 
as predictions vary in unacceptable way 1181 . 
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tiplicity distributions, which influence such ob- 
servables as rapidity distributions. It is argued 
that they should also provide us with some esti- 
mations of the fluctuations of the cross section, 
albeit we do not yet know the respective formu- 
las. However, there are already some data, which 
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Figure 5. Results on the nonmonoticity of fluctu- 
ations represented by Var(n ne g)/ (n n eg) as func- 
tion of number of participating nucleons (cf . |19| ) . 

warn us that to get such connection could be very 
difficult task. Namely, in Fig. 5 we see clearly 
that such fluctuations measured in nuclear colli- 
sions are nonmonotonic function of the number 
of participating nucleons, which in turn is pro- 
portional to the volume of interaction ^H] 8 - This 
could suggest that monotonicity observed in be- 
haviour of first moments (cf. Fig. 4) could be lost 
when going to higher moments. 
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